Low frequency response of a collectively pinned vortex manifold 
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A low frequency dynamic response of a vortex manifold in a type-II superconductor can be 
associated with thermally activated tunneling of large portions of the manifold between pairs of 
metastable states. We suggest that statistical properties of these states can be verified by using the 
same approach for the analysis of thermal fluctuations the behaviour of which is well known. The 
exponent describing the frequency dependence of a linear response is found for the generic case of 
a vortex manifold with non-dispersive elastic moduli and also for the case of thin superconducting 
film in which the compressibility modulus is always non-local. 
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I. INTRODUCTION 

In various situations the theory of weak collective pin- 
ning (for a review see Ref. treats a vortex medium 
in a type-II superconductor like an elastic manifold in- 
teracting with a random pinning potential. In particular, 
a vortex contribution to a low frequency response to an 
applied current (impedance) is ascribed to the classical 
(thermally activated) tunneling of large portions of the 
vortex manifold between different minima of the uncor- 
related random potential In the framework of this 

approach the pairs of states that allow for thermally ac- 
tivated tunneling between them are treated as current- 
biased two-level systems. Earlier analogous ideas were 
applied for the description of the dynamic response of 
spin glasses Q and of randomly pinned dislocations and 
interfaces Q. 

The frequency dependence of the two-level systems 
contribution to specific impedance Ztiiio) = Pti{'^) + 
iojlti{uj) of a type-II superconductor penetrated by ex- 
ternal magnetic field has been found by different groups 
of authors IH-Hl to be of the form: 



kiiu;) oc [ln(l/ro|c^|)]^; 
Ptiiu;)^\Lu\[\n{l/To\Lu\)]y-\ 



(la) 
(lb) 



but the conclusions of these three groups on the value of 
the exponent y are not compatible with each other. 

Comparison shows that the discrepancy appears be- 
cause the calculation of Fisher, Fisher and Huse |^ is 
based on a semi-phenomenological conjecture that a con- 
tribution of an active two-level system to inductance can 
be estimated (from above) by replacing this two-level sys- 
tem with an insulating hole of the same volume. This as- 
sumption does not take into account that the motion of 
vortices inside of a two-level system leads to a change of a 
phase distribution in superconductor (and of the energy 
of electric current in it) also outside of the limits of this 
particular two-level system. Accordingly, it turns out to 
be in direct contradiction with the explicit expression for 
a two-level system contribution to impedance derived by 
Koshelev and Vinokur and used in their calculation 



of Zti{uj), and naturally leads to a different value of y. 
The less pronounced descrepancy between the results of 
Ref. 1^ and Ref. can be explained, in particular, by 
the difference in assumptions on statistics of metastable 
states. 

In the present work we revisit this problem following 
the more reliable (as is demonstrated below) approach of 
Koshelev and Vinokur . We start by deriving the value 
of the exponent y for the generic case of a vortex manifold 
with a single elastic modulus (but with arbitrary dimen- 
sionality). Any more complex system with different (but 
local) elastic moduli will be characterized by the same 
value of y. In Ref. [Q the value of y in non-dispersive 
system has been found only for the particular case of a 
single vortex pinning. 

We then suggest a new independent way to check the 
validity of the numerous assumptions involved in calcula- 
tion of the two-level systems contribution to a linear dy- 
namic response by using the same set of assumptions for 
calculation of another quantity (the amplitude of thermal 
fluctuations), which on the other hand can be calculated 
exactly, because in the framework of a random manifold 
description it has to be the same as in absence of disorder 
. It turns out that application of the assumptions used 
in the previous calculation indeed produces an answer 
which is in agreement with the well known result for the 
pure case. 

The possibility of the additional check turns out to 
be very useful when we address the case of a thin su- 
perconducting film in which the compressibility modulus 
cii of vortex manifold is always non-local. It allows to 
draw some conclusions about size and shape distribution 
of two-level systems which otherwise would be unavail- 
able. Inclusion of these conclusions into calculation leads 
[for cii{q) oc l/q^] to very weak (double logarithmic) fre- 
quency dependence of lti{ijj) corresponding to y — Q. 

The case of a thin superconducting film is also of a 
special interest in relation with recent experimental in- 
vestigation of z{uj) in ultrathin YBa2Cu307 films 
Note that in a two-dimensional (2D) geometry with com- 
plete penetration of the magnetic field the frequency- 
dependent specific (sheet) impedance directly determines 
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the response of a film to external field, whereas in the case 
of a bulk superconductor it has to be extracted from the 
surface impedance. 

In all the cases we have considered lti{i^) diverges for 
uj 0. Thus the important consequence of our results 
is that the random manifold approximation can not be 
sufficient for the description of a truely superconducting 
vortex glass phase with finite superfluid density. 

The outline of the article can be summarized as fol- 
lows. In Sec. II the random manifold problem is briefly 
introduced and some situations when it is appropriate for 
the description of a vortex manifold in a superconductor 
are specified. In Sec. Ill the statistical properties of the 
metastable states which have to be taken into account in 
the framework of the two-level system approach are dis- 
cussed. In Sec. IV the application of this approach to cal- 
culation of the vortex medium contribution to impedance 
is presented in the systematic form for the non-dispersive 
case. 

In Sec. V the same approach is used for the analysis of 
thermal fluctuations. We show that (with the same set 
of assumptions as used earlier in calculation of a linear 
dynamic response) it indeed produces an answer which 
is consistent with expectations based on consideration of 
pure system. 

Sec. VI is devoted to the discussion of a thin supercon- 
ducting film in which the non-locality of the compress- 
ibility modulus is always important. We show that if one 
assumes that the dominant two-level systems in this case 
are strongly anisotropic (as suggested by the energy bal- 
ance estimates used in the analysis of the non-linear creep 
[pT|-Jl4[), the expression for thermal fluctuations ampli- 
tude turns out to be convergent in contrast to its loga- 
rithmic divergence in the pure system. The only way to 
resolve this contradiction consists in assuming that the 
statistics of metastable states is dominated by the pres- 
ence of hierarchical sequence of quasi-isotropic two-level 
systems. The same distribution is then used for the calcu- 
lation of the components of zti{uj). In Sec. VII the results 
are summarized and discussed. 



II. RANDOM MANIFOLD PROBLEM 

An elastic manifold (with internal dimension D) inter- 
acting with a random pinning potential can be described 
by the Hamiltonian: 

= d^Xi J d^X2[Go-l(xi-X2rV(xiK(x2) 

+ J d^x v[x, u(x)] (2) 

where the A^-dimensional vector u(x) = ii°(x) is the dis- 
placement of the manifold. The first term in Eq. ^ de- 
scribes (in a most general form) the elastic energy of the 



manifold and the second the energy of its interaction with 
a random pinning potential w(x, u). 

The simplest assumption would consist in assuming 
that the random potential w(x, u) has a Gaussian distri- 
bution with 

(z;(x,u))d = (3) 

(w(xi, Ui)w(x2, U2))d = ^(Xi - X2)w(ui - U2) (4) 

Here and further on the angular brackets with subscript 
d stand for the average over disorder, and with subscript 
th for the thermal average. We discuss only the case of a 
short-ranged random potential correlation function w(u). 

In the simplest situation (which would imply, in par- 
ticular, full isotropy and absence of dispersion) the first 
term in Eq. (|^) can be chosen in the form 

with a single elastic modulus J. 

The physical systems which can be described by 
the Hamiltonian of the form ^ include, in particu- 
lar, a domain wall in a 2D or 3D Ising type ferromag- 
net/antiferromagnet {D ~ 1,2; N — 1); a dislocation in 
a crystal (D = 1, iV = 2); a single vortex line in a large 
area Josephson junction (D = 1, A'^ = 1) or bulk super- 
conductor {D — 1, = 2); a vortex medium in super- 
conducting film {D — 2, N — 2) or bulk superconductor 
{D ~ 3, N ^ 2), a layered superconductor with in-plane 
field (13 = 3, = 1) or a large area Josephson junction 
with in-plane field (D = 2, = 1). In all these cases 
the random pinning potential is automatically provided 
by impurities present in any solid or / and by geometrical 
inhomogeneities . 

Note, however, that the random manifold approxima- 
tion assumes the energy of the interaction with the in- 
homogeneities to be uncorrelated for different displace- 
ments, whereas the energy of the interaction of an ideal 
vortex crystal with the inhomogeneities does not change 
if the vortex crystal is shifted as a whole by one lat- 
tice constant. Thus, the area of applicability of the ran- 
dom manifold approximation for the description of vortex 
crystal pinning is resticted. One can use this approach if 
the relevant displacements do not exceed the lattice pe- 
riod (Larkin regime [T5| ) or when the ordering in the vor- 
tex crystal is destroyed by the presence of defects whose 
motion with respect to the vortex manifold is dynami- 
cally frozen in comparison with the motion of the mani- 
fold itself. 

Recent experiments of the Neuchatel group on ultra- 
thin YBa2Cu307 films have demonstrated a crossover 
to the regime in which the contribution to resistivity as- 
sociated with the motion of point-like defects is negli- 
gible in comparison with the contribution which can be 
ascribed to collective pinning behaviour. 
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III. LOW ENERGY METASTABLE STATES 

The low-frequency dynamics of a weakly pinned elastic 
manifold can be associated with the thermally activated 
tunneling of large domains of the manifold between dif- 
ferent minima of the random potential In a simple 
system with a discrete spectrum only the tunneling be- 
tween the ground state and the first excited state is of 
importance at low temperatures, thus it can be reduced 
to a two-level system. For an infinite manifold one should 
take into account that such two-level systems appear at 
all scales and form a hierarchical structure, i. e. if some 
domain of the manifold can tunnel between some states 
one also has to consider the possibility of tunneling of 
smaller domains inside this area between different pairs 
of states. 

Each of such two-level systems can be characterized 
by its (linear) size L, its volume V ^ (untill speci- 
fied we discuss the simplest isotropic case with a single 
non-dispersive elastic modulus), the typical vortex dis- 
placement between the two states u, which, for example, 
can be defined by 



(6) 



the energy difference between the two states A and the 
energy barrier U which has to be overcome for moving 
the manifold (the vortex bundle) from one of the two 
states to the other. 

The universality hypothesis introduced by loffe and 
Vinokur ||] suggests that for each length scale L there 
should exist only one relevant energy scale E{L) such 
that, in particular, the typical values of A and U for a 
system of size L are proportional to E{L). The magni- 
tude of E{L) can be then estimated by estimating the 
elastic energy associated with the displacement of the 
vortex bundle of the size L: 



E(L) - JV{L) 



L2 ■ 



(7) 



If the scale dependence of the typical displacement u{L) 
is given by u{L) oc (where C is usually called the wan- 
dering exponent) Eq. (0) leads to E{L) cx with 



(8) 



On the other hand Fisher, Fisher and Huse have sug- 
gested that the scale dependence of the typical energy 
barrier U{L) can be described by another exponent il) 
not necessariliy coinsiding with x (V' ^ x)- For the sake 
of generality in the following we will keep (for a while) 
the separate notation for ■0, although we will assume that 
the typical value of A (for the given length scale L) can 
be estimated with the help of Eq. (^. 

Various properties of the manifold depend also on the 
form of the size distribution function v{L) of the two- 
level systems. The hierarchical distribution implies that 



any two-level system can include smaller two-level sys- 
tems whose size differ from that of the "parent" two- level 
system by some numerical factor of the order of one. Thus 
the ratio of the typical "neighboring" length scales has 
to be more or less constant across the whole length-scale 
range involved . This is compatible with a uniform dis- 
tribution of the logarithms of the length scales. However, 
for each length scale one should also include the factor 
1/V{L) proportional to the largest possible concentra- 
tion of non-overlapping two-level systems of size L. Thus 
v{L) has to be of the form 



dL v{L) (X 



dL 1 
~TV{L) 



(9) 



Koshelev and Vinokur Q have introduced the first fac- 
tor (imposed by the hierarchical structure) in the r.h.s. 
of Eq. (H) in the form dU/U (without any explanation), 
which for U{L) algebraically dependent on L is equiva- 
lent to dL/L, whereas in Ref. the hierarchical nature 
of the size distribution of two-level systems has not been 
taken into account. 



IV. TWO-LEVEL SYSTEMS CONTRIBUTION TO 
IMPEDANCE 



The contribution of the two-level systems to the spe- 
cific impedance of a superconductor is given [|j by 



(Vu) 



lOJ 



T \cosh^(A/2r) l + iTUj 
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(10) 



where n is the concentration of such systems, B is the 
magnetic induction and T is the temperature, whereas 
V^, u, A and t are parameters characterizing a particular 
two-level system: V is the volume in which the vortices 
are displaced (or the area in the case of a 2D supercon- 
ductor), u is the average displacement of the vortices in- 
side the bundle in the direction of current-induced force 
and A is the difference in energy between the two states. 
The relaxation time r describing the rate of the thermally 
activated (incoherent) tunneling between the two states 
depends exponentially: 

T = roexp([//T) (11) 

on the barrier U separating them. 

If one splits zti{uj) into real and imaginary parts: 



zti{uj) = pt/(^) + i^^kii^) 



(12) 



and assumes that the average over disorder can be esti- 
mated by taking for each length scale the typical (scale- 
dependent) values of all the parameters involved, the ex- 
pression for the two-level systems contribution to specific 
impedance lti{ijj) is reduced to 
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T 



dLiy{L)V\L)u'{L) {cosh- \A/2T))d{L) 

(13) 



Due to the exponentially fast increase of r with L, instead 
of including in Eq. (|lj) the factor 



1 



1 + [t{LM 



(14) 



the integration in it is cut off (at the upper limit) at a 
frequency-dependent length scale defined by the rela- 
tion t{L^)lu - 1. For U{L) oc e(L/Lc)'^ 



Lut (X Lr 



-In 4, 

e Ta\uj\ 



i/V 



(15) 



The integration interval in Eq. ( |l3| ) is limited from 
below by the (temperature dependent) collective pinning 
length Lc which determines the boundary between the 
different regimes of fluctuations. At length scales lower 
than Lc the manifold can be considered as fluctuating 
within one of the minima of the (thermally renormal- 
ized) random potential, whereas for larger scales only the 
jumps between different valleys of the potential are of 
importance. The contribution to Z(w) from length scales 
smaller than Lc has a finite limit for w — > 0. 

Since A is the difference in energy between the two 
spatially separated states in the uncorrelated random po- 
tential, the distribution function p( A) can be expected to 
remain finite for A — ^ (the broad distribution assump- 
tion d). The last factor in the r.h.s. of Eq. (|l|) (the 
fraction of "thermally active" two- level systems) can be 
then estimated [H,p| as 



(cosh-2(A/2r))d(i) 



T 



A{L) 



(16) 



where A (L) is the typical value of A for the length scale L 
[for example the width of p(A)]. Note that (for any scale) 
only small fraction of two-level systems is assumed to be 
not frozen and therefore involved in linear dynamic re- 
sponse (or thermal fluctuations). Thus they are expected 
to be well separated from each other, which justifies ne- 
glecting their interaction. 

Substitution of Eq. (|l|) into Eq. (^ leads to 



lu{uj)r^l dLv{L) 



VHL)u'{L) 
AiL) 



(17) 



According to the universality hypothesis Q A (L) has to 
be of the same order of magnitude as the elastic contri- 
bution to energy estimated in Eq. (^, which for u ^ u 
gives: 



V^{L)u^{L) 
A{L) 



V{Ly 



L2 
J 



(18) 



Substitution of Eqs. (^ and (|T|) into Eq. (|^ then leads 
to 



dL ^ 1 2 fT 1 
lti(uj) (X — —L cx-L^ — In — — 
J Jl, L J \e To\uj\ 



, (19) 



where y = 2/'0. 

With the same assumptions that have been used for the 
derivation of Eq. (|l^) the two-level system contribution 
to the resistivity is given by 



.X 1 f°° T-(L)UJ^ 



(20) 



Alternatively ptii^^) can be restored from lti{oj) with the 
help of the simplified form [|l^,0 of the Kramers-Kronig 
relation: 



2 dlnlwl 



(21) 



which is applicable for /t/(w) 
give 



/(In 1^1). Both methods 



Pt/(w) oc 



In- 



(22) 



It can be shown that y is equal to 2/'0 not only for 
the simplest case of a single elastic modulus, but for the 
general case of non-dispersive moduli. In a bulk super- 
conductor at large enough scales (which corresponds to 
low enough frequencies) all elastic moduli become local. 
The case of a thin film in which the strong dispersion of 
the compressibility modulus is unavoidable is considered 
in Sec. VI. 

For finite values of the current density j Eq. ([l^), 
which has been the starting point of our calculation, is 
applicable only if VuB/c is small in comparison with 
temperature. This determines the current dependence of 
the length scale Lj oc (T/j) -d+c at which the integra- 
tion in Eq. (|l3|) should be cut off if Lj ^ L^j- In that 
case the growth of /f/(w) with decreasing lu saturates at 
ki{uj = 0,j)x{T/j)T^. 



V. COMPARISON OF TWO APPROACHES TO 
CALCULATION OF THERMAL FLUCTUATIONS 
AMPLITUDE 

In the present work we suggest an independent way 
to check the validity of the two-level system approach 
for the description of the linear dynamic response of a 
collectively pinned manifold. This can be done because 
in the random manifold problem the static irreducible 
correlation function 

x[u\^2) - {u\^2))A)d (23) 
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(which can be associated with thermal fluctuations) ac- 
cording to Shultz et al should be exactly the same as 
in the absence of disorder: 



(K(XiK(x2)))=rGS''(xi-X2). 



(24) 



(a brief derivation can be found in Appendix). An anal- 
ogous relation for the case of periodic behaviour of w;(u) 
with respect to displacement has been suggested by Dot- 
senko and Feigel'man [T^ . 

In the presence of disorder the long-distance behaviour 
of the correlation function ( |2^ ) has to be mediated by the 
two-level systems. Therefore, the investigation of thermal 
fluctuations in terms of the two-level system approach 
and comparison of the result with well known result for 
the pure system allows to check the validity of differ- 
ent conjectures involved in the calculation of a linear dy- 
namic response. Instead of considering the dependence of 
((m''(xi)m''(x2))) on |xi — X2I one can alternatively inves- 
tigate the dependence of {{u^)) on the size of the system 
Lq, which also has to be the same as in the pure case. 

For a single two-level system with the energy gap A the 
amplitude of the thermal fluctuations of the displacement 
is given by 



{{u^))^{{u-{uUfU = 



4cosh^(A/2T) 



(25) 



For a collectively pinned manifold the dominant large- 
scale contribution to ((u^(x))) should come from the two- 
level systems which include the point x and [on the same 
assumptions as have been used while calculating Zt;(w) 
and pti{ijj)] can be estimated as 



{{u 



D) 



T / dLv{L) 



La 



A{L) 



(26) 



where the upper limit of integration is now imposed 
by the size of the system and [in accordance with Eq. 
(O)] (cosh~2(A/2r))d has been already replaced with 

t7a(l). 

Substitution of Eqs. ^ and (|) into Eq. (|2|) then gives 



dL 



L 



D-l 



(27) 



which, if compared with the trivial result for the case 
without disorder (for which ((u^)) = (u^)th): 



T 



J J\q-\>^/Lo (Stt)-" q 



d^q 1 
U3 



(28) 



reproduces all its important features. Namely, the r.h.s. 
of Eq. ( p?!) (i) contains the correct prefactor T/J; (ii) 
demonstrates the correct dependence on the size of the 
system: 



2-Z5 



L 



2-D 



2-D 



D <2 
D = 2 
D>2 



(29) 



and (iii) the system-size dependent contribution to {{ufi)) 
does not depend on the unknown disorder-related param- 
eters Lc and e. 

This allows to conclude that different assumptions 
which have been used while calculating lti{uj), pti{oj) and 
{{ufi)) [the universality hypothesis, the hierarchical dis- 
tribution of two- level systems, the broad distribution as- 
sumption for p(A)] were indeed chosen in a reasonable 
way. 



VI. THIN SUPERCONDUCTING FILM IN 
PERPENDICULAR MAGNETIC FIELD 

The long-range interaction of vortices makes the 
compressibility modulus Cn of a bulk superconductor 
strongly non-local for the wave-lengths smaller than the 
magnetic field penetration depth A. In a thin super- 
conducting film the penetration depth A is strongly in- 
creased in comparison with that of a bulk supercon- 
ductor: A = 2\^/d [jl^, where d <^ X is the thickness 
of the film. Therefore, in a thin film the dependence 
cii(<z) ~ cii/q^ (where cn « B'^/2'kA oc cgea^^) result- 
ing from a non-screened vortex- vortex interaction holds 
in a much wider range of length scales than in a bulk su- 
perconductor. Here cge ~ ^qB /32tt^A is the shear mod- 
ulus of the film (which in contrast to cn is always local), 
$0 — hc/2e is the flux quantum and — B/^q is the 
vortex density. 

If one tries to shift a vortex bundle in such a system (in 
search of the next potential minimum), it turns out that 
the optimal shape of the bundle is strongly anisotropic 
[prj-p^ with the size in the direction of the displacement 
L much larger than the size in the perpendicular direc- 
tion L± . The optimal relation between L and L± can be 
found by minimizing the total elastic energy for the given 
area of a bundle S oc LL± . Minimization (for fixed S) of 



Econ, + where 00 



— \ 2 

u 



L 



(30) 



and 



EsY. - CqqS [ — ) (31) 



or a simple comparison of -Ecom with _Esh for u ^ u gives 

(32) 



i oc ^ > i_L- 



Note, however, that when the energy scale defined by Eq. 
( ^ ) is used as an estimate for A, the factor Su^/A in 
the 2D version of Eq. (Eq) is reduced to 



(33) 



A C66 
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for arbitrary relation between L and L± . 

The calculation of Sec. V has confirmed that the form 
of the size distribution of two-level systems can be cor- 
rectly estimated by assuming that they do not overlap 
with each other, but can be situated inside each other 
forming hierarchical structures. Such estimate is consis- 
tent with the conjecture that the number of metastable 
states to which a particular domain of a manifold can 
tunnel (without moving a much larger part of the mani- 
fold) is always of the order of one Q . In what follows we 
assume that the same property holds also in presence of 
dispersion. 

The most optimistic estimate for v{L) can be then 
obtained by assuming that for all scales the strongly 
anisotropic two-level systems are arranged in the most 
advantageous way to cover all the area available, which 
corresponds to 



dL v{L) 



dL 1 

TTlZ' 



(34) 



The more realistic estimate should probably take into ac- 
count that independent anisotropic two-level systems are 
likely to have uncorrelated orientations, so the require- 
ment of non-overlapping will lead to v{L) oc L^^ . 

After substitution of Eq. ( |3^ ) into the integral [of the 
form (p6|)] defining {{u^i)) one obtains that for L oc L\ 
it is convergent at the upper limit: 



((4)) 



T 

C66 



dL 
T 



L 



< oo 



(35) 



even for the optimistic form of i^(i) given by Eq. (Q). On 
the other hand, we know that in a pure system the con- 
tribution of the transverse modes (which depends only 
on the shear modulus which is local) leads to the loga- 
rithmic divergence of Uu^)) = (u^)th- It follows from the 
results of Schulz et al 0] that in the presence of disorder 
the same behaviour should be mediated by the large-scale 
two-level systems. 

A plausible way to explain the logarithmic divergence 
of ((uj;)) consists in assuming that the film should con- 
tain not only anisotropic two-level systems with L ^ L± , 
but also an hierarchical sequence of quasi-isotropic two- 
level systems in which L± is of the same order as L. For 
such two-level systems the requirement of the balance 
between the different contributions to the elastic energy 
{Ecom ^ Esh) leads to u oc {a/L)u 4C w, which means 
that the displacement of the vortices in quasi-isotropic 
bundles is mostly of rotational type. 

The linear dynamic response has to be associated with 
the same degrees of freedom as are taken into account in 
the calculation of the static thermal fluctuations. How- 
ever, zti{uj) can not be obtained by direct application of 
the fluctuation-dissipation theorem to {{u^i)), since these 
two quantities include different linear combinations of the 
degrees of freedom involved. Nonetheless, when calculat- 
ing zu{uj) one should take into account the same set of 



two- level systems as for the calculation of ((uj;)), in con- 
trast to the case of the non-linear creep for which the 
shape of the moving vortex bundles is imposed by the 
applied current 

Application of the expression (^0|) for i?com as an es- 
timate for A shows that the factor S'^It^/A in the 2D 
version of Eq. (17) does not depend on Lj_ and can 
be estimanted as L'^ /cn ~ L^A/B^. For the hierarchi- 
cal sequence of quasi-isotropic two-level systems with 
i^iL) oc this leads to the extremely weak (double 
logarithmic) frequency dependence of 



lti{i^) oc 



7A 
52 



dL A , /T , 1 
— oc — In - In — —r 

L rJ \e To|w| 



(36) 



which can hardly be expected to be resolvable from the 
background superfluid contribution Iq in the experiments 
probing the low-frequency response of thin films. 

However, substitution of Eq. (pq) into Eq. (|2l|) gives 



Ptiii^) oc 



A 



ln(l/ro|L^|) 



(37) 



which, in contrast to Eq. (|3q), does not exhibit any spe- 
cially weak dependence on to. Note that two unknown 
disorder-related parameters Lc and e as well as the mag- 
netic field dependence have dropped out from Eq. (^7|). 

The presence of the hierarchical sequence of quasi- 
isotropic two-level systems still leaves enough place for 
more optimal anisotropic two-level systems with L± <^ 
L. Although they do not contribute much to {{u^i)), their 
contributions to the components of Zti{uj) could be of im- 
portance. However, their size distribution will be forced 
by the presence of hierarchical sequence of quasi-isotropic 
two-level systems to be of the same form v{L) oc L~^ and, 
therefore, their contribution to Iti(i^) and pti(o->) will be 
of the same form as given by Eqs. (pq)-(p^). 

In a thin superconducting film the compressibility 
modulus cii is non-local not only for Kq ^ 1 where 
cii('?) ~ cii/g2, but also for Kq <^ 1 where cii(g) w 

jl-Kq. An analogous calculation for such form oic\\{iq) 
produces for the components of Zti(io) the answers of the 
form (1) with y = l/V'- 



VII. CONCLUSION 

In the present work we argue that thermal fiuctuations 
of a collectively pinned vortex manifold are determined 
by the same degrees of freedom (related to thermally acti- 
vated tunneling between the pairs of low- lying metastable 
states - two-level systems) as its low-frequency linear dy- 
namic response. Therefore one can use the known depen- 
dence of thermal fluctuations amplitude on the size of the 
system (which has to be exactly the same as in absence 
of disorder |^|) for checking the consistency of the con- 
jectures which are used in the calculation of the linear 
dynamic response. The set of the assumptions which are 
necessary to produce the correct answer for the thermal 
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fluctuations amplitude includes, in particular, the con- 
jecture on hierarchical distribution of two-level systems 
(which means that they can be situated inside each other) 
and also the universality hypothesis. If the same set of as- 
sumptions is used for the calculation of a vortex manifold 
contribution to impedance its frequency dependence (in 
absence of dispersion) is given by Eqs. (P with y — 2/%j;^ 
where '0 is the exponent (which depends both on D and 
N) describing the scale dependence of the typical energy 
barrier U{L). 

The same result is also applicable in the limit af small 
fields, when one can neglect the interaction between dif- 
ferent vortices and treat each vortex separately as ID 
manifold. In that case one should take the value of 
corresponding to D = 1. 

Note that in the framework of our analysis the uni- 
versality hypothesis has been used only for the estimate 
of A(L). If (as suggested by loffe and Vinokur |^) it 
is further assumed that the same energy scale can be 
used for the estimate of U (L) , the value of ip will coin- 
side with X given by E q. (p[). Different approaches includ- 
ing scaling arguments |11| , |21| ] , functional renormalization 
group ]2^ , p3t and a self-consistent calculation incorporat- 
ing replica symmetry breaking [ p4| lead to 

with 1/2 < P < 1. For the case of thin superconducting 
film (D=2, N=2) or bulk superconductor (D=3, N=2) 
Eqs. (H) and (^8|) give the values of x in the interval from 
2/3 to 7/5, that is around 1. 

Koshelev and Vinokur Q have found the value of the 
exponent y only for the particular values of C and x corre- 
sponding to Z? = 1 and N = 2 and not in the general form 
as above. In Ref. |3j the hierarchical nature of the size dis- 
tribution (of the two-level systems) has not been taken 
into account and the estimate of pti{<^) has been obtained 
without integrating over the scales, which has led (in our 
notation) to y ~ 2/ ijj + \. As has been already mentioned 
in the Introduction the analogous calculation in Ref. ||] 
has been performed using the assumption which is in con- 
tradiction with Eq. ( p^ ) and therefore can not be used 
for comparison. 

Note that luiio) diverges in the limit of 0, which 
corresponds to supression of superfluid density [inversely 
proportional to \\m^^Ql{u))]. Thus the results of this 
work are not in agreement with the popular point of 
view that random manifold approach provides an ex- 
haustive description of dynamic properties of a truely 
superconducting vortex glass phase (which is supposed 
to be formed at low temperatures due to pinning | |lO| , pT|] ), 
at least if superconductivity is understood as the ability 
to carry a superconducting (non-dissipative) current and 
not only as the vanishing of the linear resistance. Our 
analysis suggests that in the framework of random man- 
ifold approach the finite value of /(/(oj — > 0) is incompat- 
ible with the correct scale dependence of {{ufi)). There- 
fore one has to conclude that the accurate description 



of a vortex glass phase which can carry a superconduct- 
ing current (if such phase exists at all) requires a more 
sophisticated treatment than the description of vortex 
medium in terms of an elastic manifold interacting with 
a random potential. For example, it possibly should take 
into account that some of the defects of a vortex lattice 
are generated by a disorder and can not freely move with 
a vortex manifold. 

However, both in the case of Larkin regime and in 
the case of dynamically frozen thermally excited defects 
(frozen vortex liquid regime) the frequency range of the 
applicability of the random manifold description of a vor- 
tex medium in a superconductor is anyway restricted 
from below. Moreover, the two-level system contribution 
to impedance lti{oj) is only logarithmic in uj and in a 
practical situation may be negligible in comparison with 
"bare" impedance Iq down to exponentially low frequen- 
cies. 

On the other hand pti{oj) produces in the low frequency 
limit the dominant contribution to the resistivity in com- 
parison with the contributions related with the normal 
channel conductance and with the oscillations of mani- 
fold within each minimum of a random potential, both 
of which at low frequencies are proportional to LlP' . 

In thin superconducting films the compressibility mod- 
ulus of vortex manifold is always non-local. This leads to 
the strong anisotropy of the vortex bundles participat- 
ing in the non-linear creep |pl]-[r^ . However our analysis 
has shown that the correct length-scale dependence of 
thermal fiuctuations amplitude requires the presence of 
a hierarchical sequence of quasi-isotropic two-level sys- 
tems. A linear dynamic response (which has to be calcu- 
lated assuming that the applied current does not change 
the properties of the system) produced by the same set 
of the two-level systems corresponds to y = 1/tp for 
cii{q) on 1/q {Aq <IC 1) and to y = [with lti{Lo) still 
diverging at w but only as a double logarithm] for 
cii{q) rx 1/q^ {Aq ^1). The contribution from the more 
optimal anisotropic two-level systems can be expected to 
be of the same form. In both regimes [cii{q) cc l/q and 
cii{q) oc 1/g^] the value of the magnetic field B drops 
out from the expression for zti{uj). 

Although in thin films the dc resistivity is always finite 
due to thermally activated motion of point-like defects of 
vortex lattice (vacancies, interstitials, dislocation pairs) 
[|l2|,^, the collective pinning behaviour has been found 
to be accessible to experimental observation [^ in the 
range of frequencies/temperatures where the activated 
contribution to resistivity is too small. 
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APPENDIX: 

In the case of a random manifold described by Eqs. 
(^-(^) it is possible to show that the irreducible corre- 
lation function of the form (^3|) remains the same as in 
absence of disorder . To this end one can express this 
correlation function through the second derivative 



of the (disorder-averaged) free energy 
F = -T ( In 



(Al) 



s?-— s^— 




(A2) 



with respect to the coefficients in the auxilary (source) 
terms added to the Hamiltonian 



H = H + s>''(xi) - s>"(x2). 



(A3) 



In order to calculate F it is convenient to shift the 
variables u(x) (over which the integration in the parti- 
tion function is performed) according to 

M°(x) => <(x) = u'^(x) -I- Gf{x - xi)s? - Gf{y: - Xa)^^ 

(A4) 

which allows to split the non-random contribution to 
Hlu"} into two terms: i/ci{<} + E{s^,s^), the first of 
which does not depend on 2 and the second 



i?(s?,4) 



GS'(xi 



X2)S?4 



(A5) 



does not depend on w°(x). 

On the other hand the distribution function of the ran- 
dom potential ^(x, u) is not affected by the shift defined 
by Eq. (|A4[). Therefore the free energy defined by Eq. 



(A2) differs from its value for original problem (that is 
for sf = S2 = 0) only by addition of the term E{si, S2), 
differentiation of which leads to Eq. (|2^). 
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